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Prediction of Clinical Irrelevance of PK Differences 
in Atorvastatin Using PK/PD Models Derived Fronn 
Literature-Based Meta-Analyses 

R Vargo^ A Adewale^, MO Behm^, J Mandema^ and T Kerbusch^ 

To support the development of a fixed-dose combination (FDC) of ezetimibe and atorvastatin for the treatment of 
dyslipidemia, bioequivalence (BE) studies were conducted across a combined dose range (10/10, 10/20, 10/40, and 
10/80 mg of ezetimibe/atorvastatin). In the BE trials, all parameters met traditional BE bounds except for atorvastatin 
peak plasma concentration (C^^^) at two intermediate doses. Literature-based metadata analysis predicted that the 
observed difference in C^^^ between an ezetimibe+atorvastatin FDC and coadministration of these agents translates 
directly into a non-clinically significant change of <1 .2% absolute difference in the percentage lowering of low-density- 
lipoprotein cholesterol. Both FDC doses were confirmed to be clinically equivalent to coadministration in the subsequent 
clinical equivalence trials. These data suggest that modeling of dose-response relationships may be useful in predicting 
clinical equivalence, lowering cost/timelines through effective powering of studies, and predicting the effectiveness of 
new dosage formulations without the need foradditional clinical efficacy trials in regulatory settings. 



To support the development of a fixed-dose combination 
(FDC) of ezetimibe and atorvastatin, bioequivalence (BE) 
studies were conducted for all combinations of the marketed 
doses: 10 mg of ezetimibe in combination with 10, 20, 40, 
and 80 mg of atorvastatin. The results of these studies dem- 
onstrated BE for all key parameters, except for atorvastatin 
peak plasma concentration {C^^^, at the intermediate doses 
(Table 1). Modeling and simulation analyses were conducted 
to understand the impact of statin-dosing regimen and formu- 
lation on low- density- lipoprotein (LDL) cholesterol (LDL-C) 
dose response and to predict the impact of the changes in 
exposure for ezetimibe+atorvastatin FDC on LDL-C lower- 
ing efficacy. 

Clinical data provide strong evidence that the clinical effi- 
cacy (LDL-C lowering) of atorvastatin on dyslipidemia is cor- 
related more with the atorvastatin dose and less with peak 
systemic drug exposure. ^"^ Reductions in atorvastatin C^ax 
are not associated with changes in LDL-C lowering efficacy.'-'''^ 
The biological processes governing the changes in LDL-C 
occur over a range of weeks to months, whereas plasma con- 
centrations of atorvastatin peak within an hour after dosing.^ 
The differences in timescale between C^^^ and pharmacody- 
namics, along with the other information described above. 



provide a rationale for total daily exposure (area under the 
concentration-time curve (AUC)) being more predictive of 
changes in LDL-C in comparison with C^^^. 

Leveraging the strong dose-response relationship for sta- 
tin efficacy on LDL-C reduction, Mandema et al.^ developed 
a dose-response model for statin LDL-C lowering as mono- 
therapy and in combination with ezetimibe. This "comparator" 
model or metadata analysis described dose responses of six 
different statins and combinations at the study-arm level. This 
work updated the dose-response model using published clini- 
cal data from 245 statin trials in 106,808 patients (Table 2). 

To predict the potential impact of the changes in exposure 
for ezetimibe+atorvastatin FDC on LDL-C lowering efficacy, 
the differences in pharmacokinetic exposure had to be trans- 
lated into LDL-C lowering. Dose-exposure models for atorvas- 
tatin AUC and C^^^were developed from the coadministration 
reference arms of the BE studies. Assuming that the observed 
exposures (AUC and C^^^ for the FDC translated into shifts 
in dose, effective dose distributions were generated as input 
for the statin dose-LDL-C model, allowing translation from 
exposure to LDL-C lowering efficacy. 

The modeling results and publicly available supporting data 
were submitted to the US Food and Drug Administration (FDA) 
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for the approval of the ezetimibe+atorvastatin FDC.^'^° Because 
the regulatory decision was primarily based on the BE data, a 
standard in providing a chnical bridge for product quality, effi- 
cacy, and safety that has been applied to other similar FDCs, 
additional clinical data were requested by the FDA. The model- 
predicted differences in LDL-C lowering for the FDC were used 
to design two clinical equivalence (CE) trials that evaluated the 
doses that did not meet atorvastatin C^^^ BE, the results of 
which were submitted to the FDA. 

RESULTS 
BE results 

Two BE studies were conducted to cover the dose range of the 
FDC. The primary hypotheses for the BE studies were that the 
observed geometric mean ratios and the 90% confidence inter- 
vals (CIs) for atorvastatin AUCg_^, atorvastatin C^.^^, uncon- 
jugated ezetimibe AUCq j^j.^, and unconjugated ezetimibe 
^max would fall within the range of 0.8-1.25. The results of the 
BE studies, performed in the fasting state, demonstrated BE 
between a single dose of the ezetimibe+atorvastatin FDC tab- 
lets and coadministration of corresponding doses of ezetimibe 
and atorvastatin as individual tablets for all of the primary 
end-point variables, except atorvastatin C^^^ at the 10/20- and 
10/40-mg doses (Table 1). Summary statistics for the trials are 
presented in Supplementary Material SI online. 

Literature-based meta-analysis 

An updated dose-LDL-C reduction model based on publicly 
available trial data was developed for statin monotherapy and 
coadministration with ezetimibe. A total of 245 of 365 trials in 
the database met the inclusion criteria for this analysis (Table 2). 
The trials represented data from 106,808 patients and evalu- 
ated all treatments across a large range of doses. A total of 1,267 
data points are in the analysis set. Each data point is the mean 
response in a group of patients at a particular time point in an 
arm of a trial plus its associated normally distributed variance 
(i.e., study-arm level data). In pooling the data, data points were 
weighted based on sample size. 

Immediate-release, once-daily administration is the stand- 
ard statin regimen for the treatment of dyslipidemia. To assess 
the impact of peak statin exposure on LDL-C reductions, trials 
testing twice-daily dosing and controUed-release formulations 
with total daily exposures (AUC) similar to those of once-daily 



Study Highlights 



WHAT IS THE CURRENT KNOWLEDGE ON THE TOPIC? 

^ Substantial clinical data provide strong evidence that the 
clinical efficacy (LDL-C lowering) of atorvastatin on dyslipi- 
demia is correlated to a lesser extent with peak systemic 
drug exposure and is correlated more with the atorvastatin 
dose. 

WHAT QUESTION DID THIS STUDY ADDRESS? 

1^ Modeling and simulation analyses were conducted to 
understand the impact of dosing regimen and formulation 
on LDL-C dose response, predict the potential impact of 
the changes in exposure for ezetimibe+atorvastatin FDC 
on LDL-C efficacy, and calculate power for CE trials. 

WHAT THIS STUDY ADDS TO OUR KNOWLEDGE 

^ The modeling and simulation analyses were able to quan- 
titatively conclude that LDL-C reduction is not highly 
correlated with peak statin exposure, which allowed for 
optimal design of the CE trials and accurate prediction 
of outcome. 

HOW THIS MIGHT CHANGE CLINICAL PHARMACOLOGY 
ANDTHERAPEUTICS 

^ This work and similar modeling analyses leveraging relevant 
competitor information can successfully predict the effec- 
tiveness of new dosage formulations and arguably may 
negate the need for dedicated clinical efficacy trials after 
near-miss BE has been found, providing faster patient access 
to therapeutic innovation. 



immediate-release formulations but with reduced peak expo- 
sures (Cj^^j,) were also included in the analysis. Fluvastatin was 
given b.i.d. in six trials (295 patients), lovastatin was given b.i.d. 
in five trials (3,562 patients), and pravastatin was given b.i.d. in 
six trials (493 patients). Fluvastatin was given as an extended- 
release formulation in 10 trials (1,558 patients) and lovastatin 
was given as an extended-release formulation in 2 trials 
(314 patients). 

The final statin dose-response model includes a different 
potency (half-maximal effective dose (ED^g)) for each statin, 
with a similar E^^^ (maximum reduction in LDL-C) and sig- 
moidicity (n) for the statins (Eq. 2). EDjq characterizes the 
dose at which 50% of the maximal response is achieved. The 



Table 1 BE results for atorvastatin and ezetimibe in healthy volunteers across the dose range for the fixed-dose combination as 
compared with coadministration of individual components 







Estimated geometric mean ratio (FDC/coadn 


ninistration) (90% confidence interval) 


Part 

(dose, mg/mg) 


N 


Atorvastatin AUCg.^ 


Atorvastatin C„,, 

max 


Unconjugated 
ezetimibe AUCg i^^, 


Unconjugated 
ezetimibe C^3^ 


10/10 


92 


0.93 (0.86,1.01) 


0.90(0.81,0.99) 


1.03 (0.99,1.07) 


1.13(1.05,1.22) 


10/20 


95 


0.89 (0.83, 0.96) 


0.77 (0.68,0.87) 


0.97 (0.90, 1 .04) 


1.00 (0.90,1.10) 


10/40 


96 


0.96 (0.94, 1 .00) 


0.81 (0.73,0.90) 


0.98 (0.93, 1 .03) 


1.03 (0.96,1.10) 


10/80 


95 


1.11 (1.06,1.17) 


1 .09 (0.99, 1 .20) 


0.97 (0.93, 1 .02) 


0.99 (0.92, 1 .06) 



AUC, area under the concentration-time curve; BE, bioequivalence; C^^^, peak plasma concentration; FDC, fixed-dose combination. 
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sigmoidicity factor captures the slope of the change from base- 
line to maximum response around ED^q. 



■N{0,(7'/N,j,] 



■Emax-Dose" 
Dose" + ED^o,. 



ycombination ^statin yezetimibe '(-^ Y'fsh 



(1) 



(2) 



(3) 



In these equations, Y^.^ is the mean LDL-C response at time 
t for the fth arm in thejth trial. Eo.^ represents the placebo 
response at time t for the jth trial. /(O^^) is the dose-response 
model for statin and ezetimibe mono- or combination therapy 
with model parameter Qj. • . j,. The parameter reflects the ran- 
dom residual error due to between- and within-subject variabil- 
ity, assumed to be normally distributed with variance O^/N^^. o is 
the SD of the percentage change from baseline, and N^^.^ is the 
sample size at time t for the /th arm in the jth trial. The relative 
influence of individual data points was weighted by treatment 
sample size during model development by influencing the vari- 
ance around that data point. In the combination (/combination^' 
/statin Statin dose-response relationship for monotherapy; 
/ezetimibe fractional change in lipids due to ezetimibe mono- 
therapy; and Y is the interaction coefficient (Eq. 3). 

The time course of LDL-C was not described because time 
did not have a significant impact on the E_^^^^ of statins, con- 
firming that a steady- state eifect had been achieved after at 
least 4 weeks of treatment. No significant difference in E_^^^ 



or sigmoidicity was found between the statins. The maximal 
LDL-C reduction for statins was found to be significantly 
dependent on baseline LDL-C (LDL.base) {P < 0.0001) and 
triglycerides (P = 0.0005), the percentage of patients with 
coronary heart disease {P = 0.0001), and whether the patients 
had acute coronary syndrome {P = 0.0004) or heterozygous 
familial hypercholesterolemia (P = 0.0003), according to the 
following model: 



itin =£max,l +£max,2 ' log(LDL.baSe/180) + 

■log(TG.base/180)+£^^_4 -CHDyo-l- 
■(ACS=yes)+£„^.5 •(HeFH=yes) 



(4) 



Ezetimibe dose response followed the same response equation 
as that for statins (Eq. 2), with the sigmoidicity factor («) fixed 
to 1. None of the covariates were found to significantly affect the 
^max ezetimibe. The covariate effects were not tested on the 
EDjg of ezetimibe because the 10-mg dose (for which most data 
are available) is at £ 

' max 

There was no significant difference in the interaction coeffi- 
cient between statins and ezetimibe among the statins. In addi- 
tion, there was also no significant impact of baseline lipids on 
the interaction coefficient. 

In the assessment of imputation methods, there was no signifi- 
cant difference in the magnitude of response between trials that 
reported the mean or median change in LDL-C. There was also 
no significant difference between trials that reported the mean 
percentage change or trials for which the percentage change was 
imputed from the mean baseline and mean on trial lipid values. 

Between-trial variances were estimated for both E and 

max 

EDjQ to quantify the impact of between-trial heterogene- 
ity on the response of statins and ezetimibe. No significant 



Table 2 Overview of trial data for literature-based meta-analysis of statins and statins coadministered with ezetimibe by treatment 



Number of 



Baseline (mg/dl) 



Data 



Treatment 


Trials 


points 


Patients 


Dose (mg/day) 


LDL-C 


HDL-C 


TG 


Age (years) 


Placebo 


133 


227 


12,236 




174(106-316) 


47 (24-70) 


1 59 (69-623) 


56 (26-77) 


Ezetimibe 


10 


45 


2,322 


10(0.25-10) 


176(159-183) 


52(43-55) 


167(142-275) 


57(52-61) 


Atorvastatin 


99 


228 


26,837 


20 (2.5-80) 


175(108-349) 


48 (24-68) 


174(101-660) 


58(26-71) 


Atorvastatin + ezetimibe 


3 


14 


352 


15(10-80) 


182(151-182) 


51 (51-58) 


168(139-179) 


59 (59-65) 


Fluvastatin 


33 


100 


5,388 


40 (2.5-60) 


200(130-296) 


49 (37-59) 


1 64 (67-443) 


56(41-76) 


Fluvastatin + ezetimibe 


1 


2 


38 


80 (80-80) 


197(197-197) 


58(58-58) 


151 (151-151) 


51 (51-51) 


Lovastatin 


26 


96 


1 0,41 1 


40(10-80) 


188(152-310) 


47 (35-53) 


174(132-265) 


56 (39-72) 


Lovastatin + ezetimibe 


1 


9 


192 


20(10-40) 


178(178-178) 


50 (50-50) 


172(172-172) 


57 (57-57) 


Pravastatin 


67 


161 


9,482 


20 (5-160) 


191 (121-314) 


48 (36-68) 


161 (59-269) 


55 (35-76) 


Pravastatin + ezetimibe 


2 


10 


254 


15(10-40) 


176(158-179) 


50 (44-50) 


177(177-270) 


57 (53-57) 


Rosuvastatin 


43 


138 


13,619 


10(1-80) 


182(114-292) 


50(31-63) 


171 (115-471) 


57 (48-64) 


Rosuvastatin + ezetimibe 


1 


1 


239 


40 (40-40) 


189(189-189) 


49 (49-49) 


186(186-186) 


63 (63-63) 


Simvastatin 


91 


190 


1 8,884 


20 (2.5-160) 


183(116-349) 


48 (24-77) 


168 (63-425) 


56 (36-75) 


Simvastatin + ezetimibe 


14 


46 


6,554 


20(10-80) 


175(134-202) 


50 (43-53) 


169(124-223) 


58 (55-66) 


Total 


245 


1,267 


106,808 




181 (106-349) 


48 (24-77) 


168 (59-660) 


57 (26-77) 



HDL-C, high-density-lipoprotein cholesterol; LDL-C, low-density-lipoprotein cholesterol; TG, triglycerides. 
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and EDcn of the statin 



between-trial heterogeneity on the £ 
and ezetimibe dose-response relationships was found, and 
the between-trial variances for these parameters were fixed 
to zero. 

The trials testing twice-daily dosing and controUed-release 
formulations were fit to assess changes in potency (EDjg), which 
would predict the impact of the blunted peak exposures {C^^^ 
on efficacy. The ED^q (based on total daily dose) of fluvastatin 
and lovastatin was significantly lower for the b.i.d. (P < 0.0001 
for both fluvastatin and lovastatin) as compared with the q.d. 
administrations (Figure 1). The ED^q values of fluvastatin and 
lovastatin were also significantly lower for the controUed-release 
formulations (P < 0.0001 for fluvastatin and P = 0.0048 for lovas- 
tatin) as compared with the normal immediate-release formula- 
tions. In fact, the ED^q values for the b.i.d. and controlled-release 
formulations were similar for both fluvastatin and lovastatin. 
For pravastatin, there was no statistically significant difference 
in potency between b.i.d. and q.d. administration regimens 
(Figure 1). Therefore, the blunted peak exposures for b.i.d. 
dosing and extended-release formulations resulted in similar 
or better efficacy as compared with the standard clinical treat- 
ment with the once-daily, immediate-release regimen for the 
same total daily dose. 



Table 3 presents the parameter estimates of the final LDL 
dose-response model. The effect of statins in terms of maximal 
change from baseline is estimated to be -76% (95% CI: -82 
to -70) in a patient population with a mean baseline LDL-C 
of 180mg/dl, a mean baseline triglyceride level of 180mg/dl, 
and no evidence of coronary heart disease, acute coronary syn- 
drome, or heterozygous familial hypercholesterolemia. The 
effect of ezetimibe in terms of maximal change from basehne is 
estimated to be -18.4% (95% CI: -19.6 to -17.2). The interac- 
tion coefficient for the coadministration of statins and ezetimibe 
was estimated to be 0.523 (Table 3). Thus, the combined effect of 
ezetimibe and statins is not simply the sum of the two effects. At 
maximal effect for statin and ezetimibe combination treatment, 
the LDL-C reduction would be predicted to be reduced by -7% 
of the sum of the two maximal monotherapy effects (Eq. 3). 

The model was found to characterize the data from each trial 
well. The characterization of the aggregate monotherapy dose- 
response relationship across the trials in a typical patient popu- 
lation for each statin and ezetimibe is shown in Figure 2. 

Prediction of LDL-C differences 

To translate the observed exposure differences between mar- 
keted atorvastatin tablets coadministered with marketed 



Twice-daily dosing 



Fluvastatin 



Lovastatin 



Pravastatin 




-20 



-40 



-60 





■ q.d. 

• b.i.d. 


\ 













Controlled release 



Fluvastatin 



20 40 60 
Dose (mg/day) 

Lovastatin 



80 




-60 



100 150 



-20 



-40 



O 

g -60 





Formulation 




■ IR 




• ER 








""^^^^-^tr —L 







-40 - 



-60 





■ IR 
• ER 







50 100 150 0 

Dose (mg/day) 



20 



40 



60 



80 



Figure 1 Impact of dosing regimen and formulation on low-density-lipoprotein cholesterol (LDL-C) dose response. A once-daily (q.d.) regimen is compared 
with twice-daily (b.i.d.) administration, and extended-release (ER) formulations are compared with immediate-release (IR) formulations. 
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ezetimibe tablets and atorvastatin in the FDC, an effective dose 
value was calculated. The calculated effective dose reflected a 
reduced dose associated with the observed reduced exposure 
in the BE trial. This reduced dose was then used to predict the 
reduced LDL-C lowering via the dose-response model. 

Two linear regression models, one describing atorvastatin 
log(AUC) as a function of log(dose) and another describing 
atorvastatin ^og{C^^^ as a function of log(dose), were devel- 
oped using the atorvastatin data from the coadministration 
arm (simultaneous dosing of each compound using the indi- 
vidual marketed tablets) of the BE studies testing the FDC, 
which included 10/10-, 10/20-, 10/40-, and 10/80-mg doses. 
The AUG and C^^^ were found to be nonlinear in the dose- 
exposure model, with slopes (95% CIs) of 1.07 (1.01, 1.13) and 
1.31 (1.24, 1.37), respectively. The slopes were then used to cal- 
culate the effective atorvastatin dose that corresponded to the 
ezetimibe+atorvastatin FDC atorvastatin exposure data for each 
atorvastatin dose. LDL-C response was then predicted from the 
effective dose distributions. 

Simulations of the LDL-C dose-response model using the 
effective dose distributions were conducted for all FDC doses 
with the three dose-exposure models: AUC, C^^^, and com- 
bined AUC and C^^^ models (Table 4). The combined AUC 
and C^^^ model is a weighted average of the resulting effective 



Table 3 Final statin LDL-C dose-response model parameter 
estimates 





Estimate 


Parameter 


(95% CI) 




-0.758 (-0.81 7 to -0.699) 


WLDL.base) 


-0.1 4 (-0.2 to -0.0796) 


f™x,3(TG.base) 


0.0506 (0.01 56 to 0.0856) 


f™x,3(CHD%) 


-0.000649 (-0.001 04 to -0.000259) 




-0.1 17 (-0.1 97 to -0.0369) 


fmax,3(HeFH) 


0.127 (0.0431 to 0.21 2) 


ED50,atorvastatin (^g/day) 


15.2 (10.5 to 22) 


ED5o,fiuvastatin(n^g/day) 


347 (214to562) 


ED50,fluvastatin(b-i-d.| 


0.645 (0.53 to 0.786) 


^f'/E'-'sO.fluvastatin 




EDsojovastatin (mQ/day) 


114(75.4to171) 


ED50,lovasta,in(b.i.d.|XR)/ED3Q,„„33,3,„ 


0.59 (0.535 to 0.652) 


ED50,pravastatin(mg/day) 


145 (947 to 221) 


EDso.rcsuvastatin (mg/day) 


4.96 (3.37 to 7.31) 


ED50,simvastatin (mg/day) 


367 (25.1 to 53.5) 


N 


0.417(0.381 to 0.457) 


fmax(ezetil^ibe) 


-0.184 (-0.1 96 to -0.172) 


ED50,ezetimibe(mg/day) 


0.228 (0.0674 to 0.769) 


Y (Interaction coefficient) 


0.523 (0.371 to 0.674) 


p (Correlation within arm by time) 


0.667 (0.614to0715) 


a (Between-subject SD) 


0.154(0.147to 0.161) 



ACS, acute coronary syndrome; CHD, coronary heart disease; CI, confidence interval; 
ED^Q, half-maximal effective dose; HeFH, heterozygous familial hypercholesterolemia; 
LDL-C, low-density-lipoprotein cholesterol;TG, triglycerides; XR, extended release. 



dose from the AUC and models, in which the weights 
are the precisions of the respective estimates. The AUC model 
predicted the smallest difference in LDL-C, and the combined 
AUC and C^^^ model falls closer to the AUC model than to the 
*-"max model. The LDL-C response curve as a function of C^^^ 
across the dose range is flat (Figure 3), driven by the observed 
nonlinearity of C^^^ and doses that are near the flat part of the 
dose-response curve. These model results predict that when 
atorvastatin C^^^ is different between two formulations (e.g., 
the largest difference of a 23% reduction for the 10/20 FDC; 
Table 1), the effect of C^^^ on the percentage of LDL-C lowering 
translates directly into a change of < 1 .2% absolute difference in 
the reduction in LDL-C (Table 4). 

Powering of CE trials 

The sample size of 406 patients for the CE trials for the 10/20- 
and 10/40-mg doses was originally calculated assuming dose- 
proportional exposure decreases in the effective dose from the 
10/20-mg BE data. The resulting effective dose from the dose 
proportionality assumption was used in conjunction with the 
published statin meta-analysis^ to obtain the prediction of dis- 
tribution of treatment difference between FDC and coadminis- 
tration of the individual components. At the sample size of 406 
patients, given the predicted distribution of treatment difference, 
the probability of success of the CE study was -80%. The sample 
size was used in the 10/40-mg protocol. Leveraging updated 
statin meta-analysis with the prediction from the simulations 
using the BE studies and dose-exposure analysis allowed for a 
more accurate estimate of the treatment difference. The sample 
sizes for the 10/20- and 10/40-mg doses were reduced to 376 
and 300, respectively, while maintaining a higher probability 
of success of 90%. 

DISCUSSION 

The large amount of literature data available for the meta-analy- 
sis resulted in a robust dose-response model for LDL-C reduc- 
tions. The model captures the data across 245 studies for both 
statin monotherapy and statins in combination with ezetimibe. 
The published data on atorvastatin show that a reduction in 
LDL-C associated with the administration of atorvastatin is 
more highly correlated with the total daily atorvastatin dose 
rather than with the measurement of peak atorvastatin expo- 

Peak exposure is blunted when the same total daily dose is 
given as a split dose twice daily or as a single extended-release 
formulation. However, the modeling analysis concluded that 
the blunted peak exposures for b.i.d. dosing and extended- 
release formulations resulted in similar or better efficacy as 
compared with that of the standard clinical treatment with 
once-daily immediate-release regimen for the same total daily 
dose. Thus, consistent with the literature data, the LDL-C 
dose-response model quantitatively predicts that a reduction 
in LDL-C associated with the administration of atorvasta- 
tin is more likely to be correlated with total daily exposure 
(AUC) rather than with the measurement of peak atorvastatin 
exposure (C^J. 
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Atorvastatin Fluvastatin Lovastatin 




E 0 50 100 150 0 50 100 150 0 50 100 150 

p 



D) Pravastatin Rosuvastatin Simvastatin 




0 50 100 150 0 50 100 150 0 50 100 150 



Dose (mg/day) 
■ IVIonotlierapy • + Ezetimibe 



Figure 2 Low-density-lipoprotein cholesterol (LDL-C) dose response of statin monotherapy and combination therapy with lOmg ezetimibe in atypical patient 
population. The solid lines are the model-predicted dose-response relationships for the difference from placebo for each statin for a typical patient population 
with a mean baseline LDLof 181, mean triglycerides of 168, and mean baseline high-density lipoprotein level of 48, and 24% of patients with coronary heart 
disease and no acute coronary syndrome or heterozygous familial hypercholesterolemia. The symbols are the mean observed values with an -95% confidence 
interval (CI) for a specific daily dose across all trialsafteradjustment for differences across the trials in covariates.TheCI for some mean observations is smaller 
than the symbol size due to the large amount of data. 



Notwithstanding the above, the most stringent ("worst 
case") evaluation of the impact of efficacy differences between 
the FDC and coadministration of individual components 
implements the assumption that C^^^ is correlated with the 
effective dose of atorvastatin. This evaluation is an extremely 
"conservative approach" because it actually runs contrary to 
the clinical evidence. These model results predict that when 
atorvastatin C^^^ is different between two formulations (e.g., 
the largest difference of a 23% reduction for the 10/20-mg 
FDC; Table 1), the effect of C^^^on the percentage of LDL-C 
lowering translates directly into a change of <1.2% absolute 
difference in the reduction in LDL-C (Table 4). This absolute 
difference based on the most stringent assumption was not 
predicted to be clinically relevant relative to the predicted 
coadministration LDL-C percentage change from baseline of 
50.3-65.1% for the 10/10- to 10/80-mg dose range. Although 
there is no standard definition of a clinically relevant change 
in LDL-C, 1.2% represents a small proportion of the -6% 
LDL-C change associated with titration to the next higher 
dose of a statin and thus would be considered not relevant in 
clinical practice. When more clinically relevant metrics of a 



Table 4 Summary of predicted difference in LDL-C lowering 
efficacy for ezetimibe+atorvastatin FDC and coadministration 
of Individual tablets (coadministration - FDC, median (95% 
confidence Interval)) 



FDC 


AUC model 


Combined AUC- 
C„,„ model 

max 


C„,„ model 

max 


10/10 


0.46 (-0.33, 0.95) 


0.52 (-0.1 4, 0.93) 


0.53 (-0.15, 1.13) 


10/20 


0.50(0.21,0.80) 


0.65 (0.34,0.95) 


1.16(0.58,1.76) 


10/40 


0.24 (0.00, 0.50) 


0.36 (0.09, 0.66) 


1.12(0.46,1.81) 


10/80 


-0.65 (-1.00,-0.29) 


-0.58 (-0.98, -0.19) 


-0.42 (-1.00, 0.14) 



AUC, area under the concentration-time curve; C^^^, peak plasma concentration; FDC, 
fixed-dose combination; LDL-C, low-density-lipoprotein cholesterol. 



change in the AUCq_^ or a combined change in both AUCg_^ 
and C^^^ exist, the models predicted a smaller absolute differ- 
ence in the LDL-C reduction of <0.65%. Thus, the impact of 
the differences in bioavailability (specifically, C^^^) between 
the FDCs and coadministration of individual atorvastatin and 
ezetimibe tablets was predicted not to result in a clinically 
relevant change in efficacy. 
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Ezetimibe 10 mg + atorvastatin 10 mg Ezetimibe 10 mg + atorvastatin 20 mg 




0.0 0.4 0.8 1.2 0.0 0.4 0.8 1.2 



CO ■ 

~ Ezetimibe 10 mg + atorvastatin 40 mg Ezetimibe 10 mg + atorvastatin 80 mg 




— I 1 1 1 1 1 1 — I 1 1 1 1 1 1 — 

0.0 0.4 0.8 1.2 0.0 0.4 0.8 1.2 



GMR of atorvastatin C^^^ (FDC/coadministration) 

Figure 3 Predicted difference in low-density-lipoprotein cholesterol (LDL-C) from coadministration for ezetimibe/atorvastatin FDC tablets as a function of 
atorvastatin geometric mean ratio (GIVIR). IVledian prediction: solid line; 95% confidence interval: shaded region; GMR of the observed bioequivalence (BE) data 
and 90% confidence interval: boxes with lines. 



These model predictions are also consistent with the labeling 
of atorvastatin, which can be given with or without food.^^ In 
the food effect study for atorvastatin, C^^^ was lowered by 25% 
and AUG was lowered by 9% when dosed with evening meals. 
However, the magnitude of the reduction in atorvastatin 
and AUG did not affect the LDL-G efficacy of atorvastatin after 
15 days of dosing.^ 

As previously noted, the LDL-G differences simulated with 
the dose-response model and BE data were leveraged to 
design the GE trials. Both cost and timelines were reduced 
by decreasing the total number of patients in the two GE tri- 
als by -17%. Both doses were found to be clinically equiva- 
lent in these GE trials (10/20:-0.2 (-1.7, 1.3), 10/40:-0.2 
(-1.9, 1.4)), % LDL-G lowering least squares mean difference 
(97.5% expanded GI)). Thus, the model-based predictions 
fully overlap with the observed lack of change in clinical end 
point. Although the model-based translation of the BE results 
was not accepted in support of the FDG filing, after review 
of the GE study results, the FDG was later approved by the 
FDA.9.10 

This work highlights the predictive power of leveraging large 
amounts of publicly available data. The modeling and simulation 
analyses were able to quantitatively conclude that LDL-G reduc- 
tion is more closely correlated with total daily exposure, rather 
than with peak exposure, with high precision and confirmation 



in the GE study results. In addition, the analyses allowed for 
optimal design of the GE trials and accurately predicted their 
outcomes. It would be expected that similar models leveraging 
relevant publicly available information can also successfully pre- 
dict the effectiveness of new dosage formulations and arguably 
may negate the need for dedicated clinical efficacy trials after 
near- miss BE study results. 

METHODS 

BE studies. Two BE studies were conducted to cover all dose combi- 
nations for ezetimibe and atorvastatin. Each dose combination was 
evaluated in an open-label, randomized, two-period crossover study 
to demonstrate the definitive BE after administration of a single dose 
of the ezetimibe/atorvastatin FDC tablet and coadministration of cor- 
responding doses of ezetimibe (Zetia) and atorvastatin (Lipitor) as 
individual tablets in healthy adult subjects. Primary end points were 
atorvastatin AUCq_^, atorvastatin C^^^, unconjugated ezetimibe 
AUCg j^jj, and unconjugated ezetimibe C^^, with the true geometric 
mean ratio of the treatments falling within 0.80 and 1.25. Additional 
study details and statistical methods are presented in Supplementary 
Material SI online. 

Literature-based meta-analysis. A database of changes in lipids, 
safety parameters, imaging outcomes, and cardiovascular outcomes 
after treatment with statins, fibrates, cholesterol absorption inhibi- 
tors, niacin, or combinations thereof reported in controlled clini- 
cal trials published in the medical literature or available from the 
FDA and the European Medicines Agency websites was developed. 
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Search terms and the inclusion and exclusion criteria are presented in 
Supplementary Material S2 online, along with the list of references 
used in the database. A total of 245 trials were identified that met the 
inclusion criteria for the current analysis. Model building and simula- 
tion were conducted in S-PLUS (version 6.2.1, Insightful, Seattle, WA) 
using the GNLS and NLME packages. 

Dose-response model. The dose response for the fractional change 
(%/100) in lipid values for statin and ezetimibe monotherapy was char- 
acterized using Eq. 1. A different placebo response was estimated for each 
time point in each trial to account for variability in the time course of 
the placebo response across the trials (nonparametric placebo model). 
Base statin and ezetimibe model structures are presented in Eq. 2, and 
the interaction for coadministration of statins and ezetimibe is presented 
in Eq. 3. 

Covariates. Time, regimen, formulation, and baseline levels of LDL-C, 
high-density-lipoprotein cholesterol, and triglycerides were evaluated for 
their impact on the maximal effect of the statin dose-response relation- 
ship. Due to the skewed distribution of the baseline LDL-C, high-density- 
lipoprotein cholesterol, and triglyceride values, the log of these values was 
used in the regression analysis. The impact of age, percentage of men, race 
(Asian vs. non-Asian), percentage of patients with diabetes, and percent- 
age of patients with coronary heart disease, acute coronary syndrome, 
metabolic syndrome, heart failure, peripheral artery disease, hyperten- 
sion, proteinuria, or heterozygous familial hypercholesterolemia on the 
dose-response relationship of statins and ezetimibe was also evaluated. 
Missing covariates were imputed on the basis of their relationship with 
other covariates. The impact of the imputation method (mean percentage 
change vs. percentage change of means) and summary statistics (mean 
vs. median) on the dose-response relationship was evaluated. 

Model selection. A stepwise addition and deletion model selection 
strategy was used, and both linear and nonlinear relationships between 
the explanatory variables and model parameters were evaluated. Model 
selection was done on the basis of a log-likelihood ratio test and SE of the 
parameter estimate at an acceptance P value of 0.0 L The SE of the param- 
eter estimate for covariate relationships was also considered. If the 95% CI 
of the parameter estimate included 0, the parameter was not included in 
the model. The clinical relevance of any relationship was also considered. 

LDL-C difference. The dose-response model was used to predict the 
mean response and uncertainty in the mean response as a function of 
dose and patient characteristics. A total of 10,000 sets of model param- 
eters were sampled from the variance matrix. For each set of model 
parameters, pharmacodynamic end points were predicted to gener- 
ate the predictive distribution as a function of treatment and patient 
characteristics. The 95% CI was taken from the 2.5th to the 97.5th 
percentiles. 

The FDC that would provide a atorvastatin AUC or C„ similar to 

^ max 

the coadministration (effective dose) was calculated from the results of 
the two BE studies. These effective dose distributions allowed simulation 
with the LDL-C dose-response model to predict the difference in LDL-C 
between the FDC and coadministration of the individual components. 
To calculate the effective dose distributions, two dose-exposure models 
were developed for atorvastatin: AUCq^ as a function of dose, and C^^^ 
as a function of dose. 



log ( AUCo-^ ) = tti log(Dose) + ttj 
log(C^ax) = [3ilog(Dose)-H|32 



(5) 



Effective dose distributions were calculated for each model, where 
DosCgjjand DosCj^^^are the effective and reference atorvastatin dose, 
respectively. IR is the individual ratio (e.g., AUCj^j,j/AUCj.j,f). 



Doscp 



( ln(IR)'l 

: Doscnef exp 

slope 



(6) 



of response. The combined i^^C-C^^ estimate of effective dose was 
derived as a weighted mean of the effective doses based on ^VC^-C^^. 
The relative contribution, or weight, of each model (AUCn-C„ ) to the 

' o ' ^0 max' 

effective dose calculation is determined by the precision of the models 
respective effective dose estimates. In other words, 



Dosepff = 



with w, = se (DoseE 



WiDose 



Eff,AUC ■ 



-WjDoseEff, 



W1 + W2 



(7) 



)and W2 - 



Doscp 



In addition, a combination of both pharmacokinetic parameters was 
considered to assess whether it provided a potentially better prediction 



). Effective dose 

distributions were simulated for all three models for aU four atorvastatin 
doses in the FDC, resulting in 12 distributions. A total of 10,000 simula- 
tions were run for each distribution. The simulations bootstrapped with 
replacement the observed individual ratios from the BE studies. The 
variability in the parameter estimates from the dose-exposure analysis 
was taken into account by using the variance-covariance matrix in the 
simulations. Each simulation resulted in the average atorvastatin dose 
for the effective dose distribution. 

For ezetimibe, the dose to provide the same exposure (effective dose) 
was calculated assuming dose-proportional pharmacokinetics. To 
estimate the potential impact of this reduced effective dose, the dose- 
response relationship for the FDC was simulated (10,000 simulations per 
dose) based on the estimated dose-response relationship of the combi- 
nation therapy with atorvastatin and ezetimibe derived from the meta- 
analysis. The simulation also included the uncertainty of the reduction 
in effective dose from the BE studies. 

Calculation of the sample size for CE trials. The original sam- 
ple size of 406 patients was chosen to provide 95% power in 
order to establish that ezetimibe+atorvastatin combination and 
ezetimibe+atorvastatin coadministration are equivalent. The sam- 
ple size calculation was based on equivalence margins of ±4%, 
assumed true difference of 1.5%, SD of within-patient difference of 
12.8%, and 85% evaluability rate using two one-sided tests,^^ each 
at a 2.5% significance level. The assumption of the true difference 
of 1.5% was based on the effective dose of atorvastatin calculated 
using a dose proportionality assumption at the 10/20-mg dose for 
both AUC and C The dose was input into the model proposed by 
Mandema eta/., "which gave a prediction of -1.5% difference. Based 
on the initial predicted distribution of difference based on the dose 
proportionality assumption, the sample size of 406 patients trans- 
lates to an -80% probability of success. 

The updated sample size incorporated the predictions from the 
updated dose-response model that accounted for the greater-than- 
dose-proportional increase in exposure. The gain in precision from the 
more elaborate modeling exercise resulted in a distinct and more accurate 
predicted distribution of differences between FDC and coadministra- 
tion ofits components for both the 10/20- and the 10/40-mg doses. This 
higher precision led to a reduced sample size, an overall savings of 17% 
across the two CE studies, and a higher probability of success (90% for 
each study). 

SUPPLEMENTARY MATERIAL is linked to the online version of the paper at 
http://www.nature.com/cpt 
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